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We study the localization transition in the integer quantum Hall effect as described by the network 
model of quantum percolation. Starting from a path integral representation of transport Green's 
functions for the network model, which employs both complex (bosonic) and Grassman (fermionic) 
fields, we map the problem of localization to the problem of diagonalizing a one-dimensional non- 
Hermitian Hamiltonian of interacting bosons and fermions. An exact solution is obtained in a re- 
stricted subspace of the Hilbert space which preserves boson-fermion supersymmetry. The physically 
relevant regime is investigated using the density matrix renormalization group (DMRG) method, 

04 . and critical behavior is found at the plateau transition. 

PACS numbers: 71.50.+t, 72.15.Rn, 73.40.Hm 

& ; I. INTRODUCTION 

Q . 

The physics of electronic transport in the presence of impurities is rich and incompletely understood. Disorder 
plays a crucial role in the integer quantum Hall effect (IQHE) as it leads to plateaus in the Hall resistance as a 
function of applied magnetic field. Transitions between plateaus are an especially interesting critical phenomenon. 
Critical behavior was predicted by Levin, Libby, and Pruisken and was observed experimentally by Wei et al. j^j 
for temperatures near the critical point at absolute zero, and by Koch et al. H for different sample widths. The 
correlation length exponent v was measured to be approximately v ss 2.3 and the dynamical exponent z was found 
to be equal to unity. 

i ^ i ■ Progress toward a theoretical understanding of the plateau transition was made with the introduction of a quantum 
tunneling network model by Chalker and Coddington Q. Subsequent numerical studies H of the Chalker-Coddington 
network model yielded values for v consistent with experiment. An argument based on the semi-classical picture 
of electron motion in the IQHE suggests that the value of v is 7/3 ||, but so far analytical derivations are still 
unconvincing. Lee |Q mapped the network model to the n — » replica limit of a one-dimensional quantum SU (2n) 
Heisenberg antiferromagnet, where the plateau transition is equivalent to a transition between two dimerized phases. 
| Lee and Wang || extended this mappingto the replica limit of an associated Hubbard model. 

■ In a continuum picture Ludwig et al. [|9| mapped the IQHE problem to a model of two-dimensional Dirac fermions 
\Q " interacting via random scalar and gauge fields. They speculated that the fixed point controlling the plateau tran- 

05 . sitions was at strong-coupling and thus not directly accessible by methods of weak-coupling perturbation theory or 
' bosonization. 

There is another field theory method besides replication which permits the averaging over quenched disorder: 
supersymmetry. McKane showed that the introduction of equal numbers of bosonic and fermionic degrees of freedom 
made the partition function independent of the particular realization of disorder. Disorder-averaged observables then 
could be expressed as the Green's functions of a field theory JhJ. Efetov used this idea to reformulate the u-model 
approach to Anderson localization [|ll| and Affleck, in a saddle-point calculation fj^] , showed that the density of 
states in a uniform magnetic field and a white noise potential was non- vanishing. Zirnbaucr fl3| applied the method 
to transport properties and derived both the supersymmetric cr-model and an antiferromagnetic supersymmetric spin- 
chain. Balcnts, Fisher and Zirnbauer ]l3 recently applied this approach to the simpler problem of the chiral metal and 
exactly solved the resulting ferromagnetic supersymmetric quantum spin-chain. Magnon excitations were identified 
H | as the expected diffusive modes. 

In contrast to the supersymmetric ferromagnet |l4|,[l5| , little is known about the antiferromagnetic supersymmetric 
spin chain. In this paper we rederive it in terms familiar to condensed matter physicists and find an exact analytic 
solution in a particular limit of the model, which opens up the possibility of studying the plateau transition pertur- 
batively. We also examine the chain in the physically relevant limit by means of the Density Matrix Renormalization 
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Group (DMRG) method and demonstrate the existence of a critical point. The correlation length exponent that we 
find is consistent with previous numerical studies of the Chalker-Coddington network model. 

II. NETWORK MODEL 

To describe the transition between quantum Hall plateaus Chalker and Coddington introduced a lattice model that 
combines the two most important features of electronic transport in a magnetic field and in the presence of disorder: 
percolation and quantum tunneling. The model consists of a network of tunneling centers which form a square lattice 
in the x-y plane; Fig. [l](a). The electron wavefunction is represented by complex amplitudes z< which are defined at 
the bond centers. In the semiclassical picture of electron transport in the IQHE [|l6| the guiding center of the electron's 
orbit propagates along a contour line of the random potential. At the saddle points, where two contours pass within a 
cyclotron radius, the electron can tunnel. Due to the presence of a perpendicular magnetic field the electron acquires 
a random Aharonov-Bohm phase between two tunneling events. In the network model this phase is represented by 
random scattering matrices (S) that relate the amplitudes of the electron wavefunction on the outgoing bonds to the 
amplitudes on the incoming bonds, Fig. [l](a). The scattering matrices are unitary and are typically chosen to be 
independently Gaussian distributed. 




FIG. 1. (a) The Chalker-Coddington network model. Amplitudes of the electron wavefunction (zj) are assigned to bonds, 
while the random scattering matrices (S) are assigned to vertices of the network. 1 and denote Hall droplets associated with 
two consecutive plateaus, (b) The Hamiltonian representation of the network model with counterpropagating chiral fermions 
as the edge states of the Hall droplets; the tunneling between the edge states is random. 

An alternative description of the network model was recently proposed by Lee |?J where the random scattering 
matrices are replaced by a Hamiltonian which describes a one-dimensional array of counterpropagating edge states in 
the presence of disorder, as shown in Fig. |l](b). The direction of propagation will be later interpreted as the imaginary 
time; to emphasize this point of view we relabel the y-coordinate as r and also denote the x-coordinate with the 
integer index n. 

Expressed in terms of electron creation xjfi and annihilation ip operators the Hamiltonian may be written as 

H[4>\ i>] = ^2 [( _1 )™^,r( i5 r + uy r )?/V T ~ (*n,^i, T i+l,r + Ki,Ti>n+l,T^n,r)], (1) 

71, T 

where the sum is over the discrete vertices of the network model, and i/j^ T (id T )'ip n , T is defined in momentum (Fourier) 
space, in the limit of small k T , as ^r n k (k r )tp n ^k T ■ To avoid ordering ambiguities, for now we keep r discrete and refrain 
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from taking the continuum limit until later, after the average over disorder has been performed. The Hamiltonian, 
Eq. [I], reflects the chiral nature and linear dispersion of the edge states, with alternating propagation forward and 
backward in the r direction at Fermi velocities ±1. It also incorporates complex tunneling amplitudes t n _ T with 
random phases, which represent the random Aharonov-Bohm phases of the electrons discussed above, as well as 
random on-site potentials w n ^ T . The random terms are taken to be independently Gaussian distributed with zero 
mean, and variance 

<C W n<T W n ',T' > = 2U 6 n<n i S t ,t' 

« C,r **',T' » = Jn <W S T ,r' , J„ = J X (1 + (-l) n fl) . (2) 

Staggered modulation in the parameter J n accounts for differences in tunneling between counterpropagating electrons 
in adjacent plaquette-columns as depicted in Fig. 0(b). As R — > — 1 electrons only circulate clockwise around the 1 
plaquettes. We can picture these plaquettes as isolated lakes of electrons (regions of attractive potential with filling 
fraction 1) separated by dry land (repulsive potential and filling fraction 0). In the opposite limit R — ► 1 electrons 
circulate counterclockwise only around plaquettes, which corresponds to isolated islands (empty regions of filling 
fraction 0) in the midst of an ocean of (filling fraction 1) electrons. Precisely at R = the land and sea cover equal 
areas and we expect critical quantum percolation associated with the transition between the and 1 plateau. The 
percolating edge states are characterized by a diverging length scale, £ ~ R~ v . Critical behavior in accord with this 
scenario was found from the disorder averaged transport Green's function of the Hamiltonian (Eq. |l|) by Wang and 
Lee [||. From a Monte-Carlo treatment of the replicated Hamiltonian they obtained a correlation length exponent 
v = 2.33 ± 0.03, in good agreement with experiments and other numerical simulations of the network model. 



III. TRANSPORT GREEN'S FUNCTIONS 



To study the localization-delocalization transition in the network model as described by the random Hamiltonian, 
Eq. |l|, we focus on the disorder-averaged transport Green's function [l7| 

A"(1,2)=<G + (1,2)G_(2,1)» (3) 

where 

G±(l,2) = (ll 1 [2) (4) 
H — hi ± lrj 

are the advanced and retarded Green's functions; E is the energy and rj a positive infinitesimal. We introduce the 
usual shorthand for coordinates with 1 = (ni,7i) and 2 = (n 2 ,r 2 ). In the localized regime, R is non-zero, and we 
expect the transport Green's function to be short-ranged 

X(l,2)~exp(-r 12 /0 (5) 

where m is the distance between points 1 and 2, and £ <~ R~ v . The correlation length £ reflects the finite extent of 
the electron wavefunction. To calculate the matrix elements of the resolvent of H appearing in the expression for the 
advanced and retarded Green's functions, Eq. we introduce a pair of complex scalar fields 

«-'>-(i&3)- < 6 » 

The Green's functions jlTj] can now be written as 

G+(l,2) = -i (<Mi)<OT> 

G_(l,2) = +i <0 i (l)0I(2)> , (7) 
where averages of observables O are taken with respect to a Gaussian measure, 

(O) =Z- X J D0D0* 0[<t>,4>*} cxp(-S[0,</>*]) 
S[4>,4>*] = +if?$>T(ra,r)0|(n,r)} - ir?^^K^)€( n ' r )} ■ ( 8 ) 
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In the above equation we have set E = 0; in the advanced and retarded Green's functions non-zero energy is easily 
recovered by replacing 77 with r/±iE. The choice of signs in the terms proportional to 77 is dictated by the requirement 
that the Gaussian integrals be convergent for 77 > 0. The quadratic form H[<j>] is obtained simply by replacing the 
operators ip n ,r and T in the Hamiltonian of the network model, Eq. [I], with complex scalar fields <fi(n, r) and 
4>* (n, r) . Finally, the normalization factor 



Z = J D^)D0*exp(-S'[0,0*]) (9) 
can be written as Gaussian integral over Grassman fields [lS|] 

Z- 1 = I Di/>DV>exp(-S'[i/>,T/>]) (10) 



where the anticommuting fields 

are the supersymmetric partners of the scalar fields. Using both Eqs. [lO] and |§|, the correlation functions, Eq. (7), may 
be written as a combined Gaussian integral over scalar and Grassman fields 

(O) = J DcpD<p*D^Dip 0[<p,<p*} exp(-S[cf>,cf>*} - S[ip,if>}). (12) 

The utility of this representation lies in the fact that, for each realization of disorder, the normalization factor has 
been absorbed into the action. Now the average over disorder can be interchanged with integration over the scalar 
and Grassman fields, and the problem of computing the disorder-averaged transport Green's function, Eq. ||, reduces 
to the problem of calculating correlation functions in a pure field theory, albeit one with interacting fermionic and 
bosonic degrees of freedom. Our strategy in what follows will be to perform the disorder average in the path integral, 
Eq. [l^. Then, because the r-derivative term is linear, the disorder-averaged path integral will be interpreted as 
a coherent state path integral for a one- dimensional Hamiltonian. The problem of calculating disorder-averaged 
correlation functions for the network model which is infinite in the r direction is thus transformed to the problem of 
computing expectation values in the ground state of this Hamiltonian. 



IV. NON-HERMITIAN HAMILTONIAN 



The supersymmetry method, as discussed in the preceding section, allows us to perform the average over disorder 
inside the path integral. Before proceeding with this step we first transform the scalar and Grassman fields on even 
sites, 



n even: 



and on odd sites 



n odd: 



<P*(n,T) -> <Aj(ti,t) i>i{n,T) -> ipi(n,r) 

T (n, t) -> 0f (n, t) ?/) T (n, r) -> -$ T (n, r) 
<t>] (n, t) -> T (n, t) (n, r) Vt ( n i T ) 5 



(13) 



(14) 



the remaining fields are left unchanged. The transformation is guided by two requirements. First, the linear derivative 
term in the action acquires the canonical form for both scalar and Grassman fields |18|, 



ip(n, t) ■ d r ip(n, t) + <p*(n, r) • d T 4>(n, r) 



(15) 



Second, the sign of the n term in the bosonic part of the action must be such to ensure convergence of the path 
integral; see Eq. [lj below. 
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Next we perform the disorder average inside the path integral using Eq. ^: 

«exp(-S[</>,</>*]-S[-</>, </>])» = cxp(~5 SU sy[0,0*,-0,^]) . (16) 
The action Ssusy can be written as a sum of three terms 

Ssusy = S T + Su + Sj + S n (17) 

where S T is given by Eq. [l^, and 

Su = U ^2{4>(n,T)a z ip(n,T) + 4>*(n,T)a z 4>(n,T)} 2 

n,T 

Sj =J2 J ™ [(- 1 )" +1 C(n, t) + B(n, r)][(-l)"C(n, r) + B*{n, r)] 

n,r 

S n =V ^[V'(n,r)--0(n,T) + 0*(n,T)-0(n,T)] . (18) 

n,T 

a z — diag(l, —1) is the 2x2 Pauli matrix operating in spin space (f , J.), and we have introduced new fields bilinear 
in the original Grassman and scalar fields which live on the bonds, 

B{n,r) = (j>*{n,T)<j>*{n + l,T) - <p L (n, r)<j) L (n + 1, r) 

C(n,r) = ^ T (n,r)^ T (n + l,r) + V>].(n,T)V>j.(n + 1, t) . (19) 

Now it is safe to take the r-continuum limit and interpret the r direction as imaginary time. The action S'susy 
is precisely that which appears in a coherent state path integral representation of the partition function for a one- 
dimensional Hamiltonian given by: 

ffsusY = H v + Hj + H v . (20) 

The different terms in the Hamiltonian can be read off from the corresponding terms in S'susy only after the scalar 
and Grassman fields have been properly ordered: the fields are commuted and anticommuted so that (jf and ip 
always appear before </> and ip- The normal ordering is essential jl8| because the scalar and Grassman fields appear 
in the coherent state path integral as left and right eigenvalues of creation and annihilation operators, respectively. 
Introducing bosonic operators which satisfy the usual canonical commutation relations [fef, Sj] = [hi, St] = 1 and 

fermionic operators which satisfy canonical anticommutation relations {£-]■, ct} = {cj, cj} = 1, normal ordering, and 
dropping the hat over the operators for clarity, we obtain 

n 

Hj = ]T J n [: B n B\ : - : C n C n : + (-l)" +1 (C„^ - & n B n )] 

n 

H, = V Win + A/?„ + A/-f„ + A/|J ■ (21) 

n 

We have introduced the boson and fermion number operators A/"' = btb^, A/J = St&j, J\f^ [ = ctcj, and A/J = c\c[- The 
remaining bilinear operators take their form from Eq. |l^: 

B n = &|„&}„ +1 ^ b ln b ln+1 

C n C |n C T«+l C i" C l"+l ' (22) 

and products of these operators which appear in Hj are normal ordered. The terms Hjj and H n are Hermitian while 
Hj is not. The lack of Hermiticity has its origin in our treatment of one of the spatial directions as the time direction. 
Evolution forward in real time must be unitary, but no such restriction holds for evolution in a spatial direction. 
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V. EXACTLY SOLVABLE LIMIT 



The ground state of -ffsuSY holds the key to understanding the localization transition in the network model. The 
disorder-averaged transport Green's function can be written as a ground state average, 

K(l,2) = (4(l)4(2) C J(2)cJ(l)) SUS Y (23) 

if 1 is an odd lattice site and 2 is an even site. This equation follows from Eqs. ^ and [?], once the transformations 
in Eqs. ^ and ^ have been taken into account, and after the Grassman fields have been replaced with the fermion 
creation and annihilation operators in the Heisenberg representation. To find the ground state we first perform a 
canonical transformation 



cj. -> c{ bi -> b l = b[ 

4 -> c l ~ b l 



(24) 



on every site of the chain. This transformation preserves the canonical commutation and anticommutation relations 
between creation and annihilation operators, and it ensures that ifsusY commutes with the number operator Af n = 
Af+ n +A/ r f„ +A/?„ +A/"f„ on each site, where the new A/J = Note, however, that in the bosonic sector the new 6j_ 
and bi are not Hermitian conjugates of each other, rather, (frj.)' = — b±. The transformation also removes the explicit 
SU(2) spin symmetry as the |-spin operators are unchanged. 

As [-ffsuSY, Afn] = the Hilbert space breaks up into sectors with different numbers of particles on each site, and 
we are faced with the task of identifying the lowest energy sector. After performing the canonical transformation, 
Eq. |24|, the new Hjj term is 

Hu = U ]T[(A/- T fc „ + M* n + Aff n + A/fJ 2 - (A/?,, - N\ n + Aff n Af?J} . (25) 

n 

For large U the ground state lies in the subspace of the total Hilbert space which has either no particles, or just one 
"f-fermion or one |-boson per site, as there is a relative energy cost of 2U for a |-particle. The energy of the vacuum 
state, however, is raised with respect to the singly-occupied states once the other two terms in the Hamiltonian, Eq. ||l| 
are included. As the ground state occupancy is thus non-zero, the density of states (DOS) is also non-zero. This is as 
expected because the DOS does not vanish either at the Anderson transition (l7) or in the Landau level problem |jl2| . 
For large U the supersymmetric exchange Hamiltonian (Hj) can be written as 

Hj = Jn WAi - A/?„A/-f n+1 + (-l)" +1 ( C t T „6 T „cj„ +1 6 Tn+1 + b\ n c 7n b\ n+1 c u+1 )] ; (26) 

n 

all the terms with J,-spin fermions or bosons have been eliminated. Physically this simply means that there is no 
transport in the U J limit. 

The large-£7 Hamiltonian, Eq. |2^, is non-Hermitian, degenerate, and it is also defective: its eigenstates do not span 
the Hilbert space. The defect in Hj is evident even for just two sites with a four-dimensional Hilbert space spanned 
by: {fcj^lO), fejjC^lO), cjj&j^O), cj^cjjjIO)}. Only three (not four) distinct eigenstates of Hj can be constructed. 
Dropping the spin index for clarity, and setting J n = 1, the two site Hamiltonian may be written: 

Hj = blb 2 b\bi — c\c2c\c\ — c\c\b\b2 + b\b\cic 2 ■ (27) 

As the Hamiltonian is non-symmetric, there are both right eigenstates: 

Hj |**(a)> = e a \* R (a)) (28) 

and left eigenstates: 

(V L (a)\Hj = (V L (a)\e a , (29) 

with fi> L (a)\ = (| , I' L (a)))t ^ {^> R (a)\ in general. In this case the three eigenstates are degenerate, with zero energy. 
The right eigenstates are: 

\* R (l))=c{bl\0) 
\* R (2))=b\ci\0) 

|**(3)> = (blbl + c\4)\0) (30) 
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as may be easily verified by direct substitution into Eq. ^8|. The defect may be removed with the addition of a small, 
supersymmetry breaking, chemical potential for the fermions: 

Hj Hj + fi c N c ; /i c <0 (31) 

where N c = 5Z n (.A/+ n +-^?n) * s * ne total fermion number operator. This term lifts some of the degeneracy and splits 
the third state into two distinct ones: 

|¥ fl (3±)) = -±= [(-/x c ± y/nl - 2^)6$ + (/i c ± Vm? - 2^)44] |0) (32) 

V JV 



with energies e^± = \i c ± v/juf — 2/i c ; iV is the normalization constant. The leading dependence of the energies on 
y 7 /^, as opposed to linear dependence on /i c , is a telltale sign of the defectiveness of the Hamiltonian in the /j, c — > 
limit. The left eigenstates are easily found: 

<f L (l)| = <0|6 2Cl 
<f L (2)| = (0|c2&i 

(* x (3±)| = <0| [(- Mc ± y/tf - 2 /i c )6 2 &a - ( Mc ± - 2^)c 2 cx] ^ (33) 

Now the four states can be properly normalized, 

(^(6)1^(0)) = 6 ab , (34) 



by setting N = -4/x cV Vc _ 2/v 

Before turning to the problem of more than two sites, we examine the eigenstates of -ffsusY with |-spin excitations 
in the physical, finite-C/ limit. By direct examination of the Hamiltonian we find that eight right-eigenstates can be 
formed by placing one |-spin and one |-spin particle on the two sites: 

& + T1 4 2 |0), ^lO); (35) 

the other six are obtained by exchanging b with c or | with J,. Two more eigenstates with two |-spin excitations are 

cj^lO), bj lC | 2 |0). (36) 

These ten states, as well as two other states which are not as simple to write down, are separated by O(J) + 0(U) 
gaps from the remaining four zero-energy ground states. The existence of a gap for spin excitations means that the 
transport Green's function, Eq. S, decays exponentially in the r direction. Physically this may be interpreted [l9| ] 
as the expected one-dimensional Anderson localization which is driven by backscattering among the two counter- 
propagating edge states. 

Surprisingly, the problem of more than two sites can be solved exactly in the U 3> J limit. To see this we consider 
the |-Hilbert subspace with just two states per site which may be represented as 

&}|0)-|0) 

c{|0>-/t|0>, (37) 

where the /'s correspond to spinless fermions. States with more than one fermion per site are forbidden by the Pauli 
exclusion principle; those with more than one boson per site are disfavored energetically by H n and the first term 
in Hj, Eq. The elimination of the (zero in energy) vacuum state is justified a posteriori as the ground state 
energy will turn out to be negative; see Eq. ^2] below. The different operators appearing in the Hamiltonian (Eq. |26| ) 
projected onto the subspace spanned by |0) and /^|0) can be represented solely in terms of the new fermion operators, 

b\h - 1 - /V b\ C] - / 

(38) 

4c T - pf c|6 T -> /t 

This transformation is the analog of the Jordan- Wigner transformation |lq] which maps the x-y spin chain into a 
system of non-interacting tight-binding fermions. In this representation the Hamiltonian is quadratic, and setting 
J = 1 we obtain: 
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Hj = ~J2(l+ i-l) n R) [flfn + fl+Jn+l + (-mtifl+1 + fnfn+l)} 



(39) 



The Hamiltonian is now easily diagonalized. Due to the alternating term in Eq. |3^ it is convenient to introduce 
separate fermion operators for even and odd sites, 



fin — d n 



(40) 



where on the odd sites we have performed yet another particle-hole transformation. With this substitution, and upon 
Fourier transforming to fc-space, the Hamiltonian may be rewritten as: 



Hj = J2 [ 2d \dk - 2e{e k + (1 - R){d\e k - e\d k ) + (1 + R){e lk d\e k 



-ik 



Ad k )} 



(41) 



Diagonalization is now a simple matter of computing the eigenvalues of a 2 x 2 non-Hermitian matrix, and we obtain: 

e k = ±2y/l-R 2 | sin(/c/2)| . (42) 

The eigenvectors are easily found. The k = piece of the Hamiltonian is defective, and once again the defect is lifted 
with the introduction of small negative \x c . 

The ground state is constructed by filling up all of the negative energy states with the spinless fermions in the 
reduced Brillouin zone < k < ir. Two things are remarkable about the dispersion relation: First, no gap opens 
up for R 7^ 0. The expected spin gap at non-zero R cannot be recovered in the U 3> J limit because none of the 
excitations in the reduced Hilbert space carry spin. Second, the dispersion at the left and right Fermi points is linear, 
unlike the quadratic dispersion found for the ferromagnetic supersymmetric chain which describes the chiral metal 
fl4j| . The linear dispersion is important to recover here because in the original plateau transition problem, unlike the 
chiral metal, the two spatial directions are equivalent. 

In Hcrmitian supersymmetric theories the ground state is unique and has zero energy. Excited states pair into 
sectors of even and odd numbers of fermions and consequently the partition function Z = STr e~ l3HsvSY — 1 at 
all temperatures 1//3. Here the supertrace operator STr is simply the ordinary trace over states with the insertion 
of an additional operator (— 1) °, where N c is the total fermion number operator introduced above. The operator 
(— 1) compensates for the anti-periodic boundary condition imposed on Grassman fields in the r-direction of the 
corresponding path-integral Q . 

As the ground state energy is negative for the non-Hermitian Hamiltonian considered here, it is important to 
clarify how supersymmetry is manifest. In particular, at non-zero temperature (corresponding in the original network 
model to finite spatial extent and periodic boundary conditions in the y-direction) it is important to check whether the 
partition function really equals unity, Z = 1. We first note that each term in the Hamiltonian i/suSY = Hu + Hj + H n 
either does not change the total fermion number, or changes it by an even integer. Therefore the many-body Hilbert 
space splits into two separate sectors: one with an odd number of fermions and (— l) Nc = —1, the other with an even 
number and (— 1) = 1. The map between these two sectors is realized by fermion-valued supersymmetric charges, 



= £ b\ n c ]n - (-l) n c}„6 T , 



E 



(-l)™6 t Tn c Tll + c}„6 T , 



(43) 



which are the analogues of those found for the supersymmetric ferromagnet p4| . 

Straightforward algebra confirms that [Hsusy, Qil^r [^SUSY, Q2] = in the /i c — * limit, and also that {Qi, Q2} — 
0. This holds true for -Hsusy in the finite-?/ (Eq. |2l| ) and in the infinite-?/ (Eq. |26| ) limits. Therefore the charges 
Qi and Qi map one eigenstate with an odd number of fermions into a second state, degenerate with the first, but 
containing an even number of fermions. For example, two of the four degenerate eigen-states of the two-site problem, 
Eq. [SO], have a single fermion while the other two states, Eq. |32|, have a mixture of zero and two fermions. In this 
case we find: Qi\^ R (l)) = Q 1 \ 1 ^ R (2)) cx (|^ fl (3+)) - \^ R (3^)}). Operating with the charges a second time, on 
states ^^(Si)), then returns states | X 1/' R (1)) and \^> R (2)). In this way states which differ in content by one fermion 
are paired and also grouped into degenerate quartets. The only exception is the unique, zero-energy, empty vacuum 
state, |0), which is annihilated by both supersymmetric charges, Qi|0) = Q2IO) = 0. All the other states (with either 
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positive or negative^] energy) pair off and Z = 1 at any temperature as expected from the supersymmetry of the 
construction. For the remainder of this paper, however, we consider only zero-temperature properties. 

VI. DMRG CALCULATION 

To make further progress in the physically relevant finite-£/ case we employ the Density Matrix Renormalization 
Group (DMRG) introduced by White (20]. We utilize the simpler and faster "infinite size" algorithm. It permits 
accurate evaluation of observables in the thermodynamic limit of large chain length but is less accurate for finite size 
chains. 

As the Hamiltonian is non-Hermitian, both right and left eigenstates must be found. We employ the Davidson-Liu 
algorithm Q| extended by Morgan |22] to non-symmetric matrices. Sixteen on-site operators with a = 1, 2, • • • , 16, 
are introduced. In normal-ordered form these are: 



(44) 



After the canonical transformation (Eq. £4|) the exchange part of the supersymmetric Hamiltonian, Eq. |2l], can be 
represented in terms of these on-site operators as: 



S 1 : 


= &}6 T 


5 3 E 


= b\h 


s 9 = 




S 13 = b\c, 




= Mt 




= kh + 1 


5 15 e 


E h c 1 


s 11 = l iCl 


glO 


= 4& T 


S 16 


= 4 6 i 


S 5 = 


4 C T 


S 7 = cjcj. 


S 14 


= 4 & T 


S 12 


^4 b i 


s s = 


c}c T 


S 6 = c\ Cl - 1 



Hj — - Jn 



8 16 
a—1 a— 9 



(45) 



with coefficients A a = —1 for a = 1, 2, 10, 11, 14, 15 and otherwise X a = +1. The second term in the above equation 
is non-symmetric; it is composed of the eight operators (S a , a = 9, • • • , 16) which have fermionic character. When 
represented in terms of matrices with real- valued entries, the fermionic nature of these eight on-site operators leads 
to additional (— ) signs, which have been taken into account in writing Eq. In particular, these (— ) signs are 
dictated by the fermion ordering convention along the chain which we adopt to define the many-body basis states: 
cj • • • ct • • • c' k ■ • • |0) with i < j < k. As the Hamiltonian commutes with the number operator on each site, Af n = 

J\ff n + J\f° n + + A/J„ , the on-site operators 5° may be represented as 4 x 4 matrices acting on the 4-dimensional 
single-site Hilbert space 

|1> = 6+10) |3)=4|0) 

(46) 

|2) = 5;|0) |4)=cJ|0). 

For example, the only non-zero matrix element of S 1 is { 1 1 iS' 1 1 1 ) = 1. Note that the state (2| = (0|(6|)t = — (0|6| 
but the (— ) sign may be dropped as it always appears in pairs in the exchange Hamiltonian Eq. ^5|. We note that 
down spins are energetically disfavored due to the ±1 terms appearing in the definition of operators S 2 and S 6 . The 
constant terms originate in the canonical transformation, Eq. 

As in the analytically tractable U 3> J limit, the ground state remains four- fold degenerate in the fi c — > limit even 
down to U = 0. This is accounted for by the existence of two supersymmetric charges, Eq. [l3|, that commute with 
the Hamiltonian. The ground state expectation values of observables O, which correspond to thermal expectations 
in the (3 — > oo limit, are computed using the supertrace formula, 

4 

(0)=^ L (a)\(-l) N "0\^ R (a)) . (47) 



1 Note that H ^ Q'Q as it would for a Hermitian supersymmetric field theory; hence, negative energy states are permitted. 
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That both left and right eigenstates appear in Eq. 47 has its origin in the "resolution of identity" which is inserted 
into the partition function for the purpose of computing expectation values: 



1 = £>*(<*)) (* L (a)| 



(48) 



here the sum runs over the entire Hilbert space. (The states span the Hilbert space because the defect was eliminated 
via the introduction of p, c < 0.) 

We introduce reduced density matrices for both of the augmented blocks, each of Hilbert space size AM, by 
computing partial traces over states on the right (or left) halves of the chain. The reduced density matrix for the left 
half of the chain is given by: 



am r 



Pij 



a=l 



9§,(a) *&(a) , *&(o) *£,(a) 



(¥*(o)|tf*(o)> (* i (a)|* i (a)) 



(49) 



a similar formula defines the density matrix for the right half of the chain. Here *&fj} R = (i,i'\^ L ' R ) are the matrix 
elements of the many-body wavefunction projected onto a basis of states labeled by unprimed roman index i which 
covers the left half of the chain and primed index i' which covers the right half. The states are normalized in Eq. |^ 
to insure that Trp = J2i Pa — 1 ■ As can be verified numerically, all of the eigenvalues of p are real (as it is a 
symmetric matrix) and, more significantly, they are positive and hence may be interpreted as probabilities. The 
DMRG truncation procedure thus corresponds to the elimination of the least-important sectors of the Hilbert space. 
Our choice for the reduced density matrix agrees with that of Bursill, Xiang, and Gehring |^3j but differs from that 
of Nishino Q who argues that a product of left with right eigenstates is optimal in a variational sense. However, the 
resulting reduced density matrix is non-Hermitian and some of its eigenvalues are negative; they cannot be interpreted 
as probabilities. Further, it is a difficult numerical problem to find a large fraction of the eigenstates of a non-Hermitian 
matrix. Since the reduced density matrix is introduced solely for the purpose of truncating, in a systematic manner, 
the rapid growth of the Hilbert space with chain length, here we adopt the symmetric choice (Eq. ^) and check that 
observables do not change significantly when the size M of the block's Hilbert space is enlarged. 

The exactly solvable U ^ J limit permits another important check on the quality of the DMRG algorithm. We have 
verified that the DMRG procedure accurately reproduces the energies of the four degenerate ground states obtained 
from the dispersion relation (Eq. |42[ ) as the chain grows. For example, in the case U — 10 J and for a block Hilbert 
space size of M — 80, the energies agree to within 1 part in 10 4 at chain lengths up through L = 16. 

We continue to lift the degeneracy, and the defectiveness, of the Hamiltonian via a small chemical potential for 
fcrmions, Eq. Bl]. Typically we set \l c = — 1.0 x 10~ 6 in units where J = 1. Although each of the four terms in 
Eq. ^| in general depends strongly on \x c we have verified that this dependence drops out in the sum over the four 
nearly degenerate ground states (see Fig. ||). The rest of the DMRG algorithm is standard: During each iteration, 
the chain is cut into left and right pieces, the reduced density matrix for each half is formed and diagonalized, new 
representations of the operators are found, and finally two sites arc added to the middle of the chain and the resulting 
new Hamiltonian is diagonalized via the modified Davidson-Liu algorithm. The near degeneracy and defectiveness of 
Hj slows down the diagonalization considerably and can be partly compensated by increasing the size of the Davidson 
subspace. We run the double-precision, vectorized, multi-processor C-code on a Cray EL-98 computer. The largest 
computations take about 1 Gbyte of memory and require several weeks to run. 



VII. DMRG RESULTS 



Previous work on the replica limit of the Heisenberg [[7j and Hubbard |8| models focused on the scaling of the spin 
gap A S (R) as a function of R, the staggered coupling between neighboring sites. It is difficult to gain direct access 
to the spin gap in the supersymmetric chain due to the plethora of low-lying excitations in the pure spin-up sector 
(discussed in Sec. [yj) which persist at finite-C/. To isolate spin-flip excitations would require finding many low-lying 
eigenvectors, a problem for which the DMRG method is not well suited p3[ . 

Critical behavior is evident, however, in the expectation value of a spin-flip operator which lives on the bonds 
between neighboring sites. It is important to note though that according to Eq. |4?] supersymmetry-invariant operators 
(those which commute with the supersymmetric charges Qi^), have zero expectation value. For example, (-ffsusy) — 
despite the fact that the ground state energy is negative. The non-supersymmetric fermion spin-flip operator 
(— 1) °S„iS„ +1 , on the other hand, is an appropriate choice. Even with no applied dimerization (R = 0), its ground- 
state expectation value 
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((-1)"" S 7 n Sj l+1 ) = ((-1)^ c\ n c [n c\ n+1 c ln+1 ) ee B(n, n + 1) (50) 

has even/odd modulation due to the free boundary condition at the chain ends. This modulation can be isolated by 
defining 

A(n) ee (~l) n [B{n, n+l)-B(n- 1, n)] + c.c. , (51) 

and it can be understood as follows: chains with an even number of sites and free boundary conditions at the chain 
ends have no reflection symmetry about any site center. The chain ends induce slight dimerization (which also breaks 
reflection symmetry about site centers). Each time the chain length is increased by two sites during a DMRG iteration 
the pattern of weak and strong bonds alternates; the term (—1)™ appearing in Eq. |si] accounts for this alternation. 

As A(n) is the expectation value of an operator which induces dimerization at R =/= it has scalin g di mension xa 
which is related to the correlation length exponent by the standard scaling relation: v = (2 — xa)" 1 [p5[ . Therefore, 
in a finite system we expect 

A(n)~L" XA , L^oo. (52) 

This observation provides us with a means of extracting xa from the DMRG by measuring the expectation value of 
the dimerization operator on the central bonds, A(L/2), at each DMRG step. 

To test this method, and the DMRG algorithm, we first examine the ordinary isotropic, nearest-neighbor, spin- 1/2 
SU(2) quantum antiferromagnet. In this case the dimension of the dimerization operator (—1)™ S n ■ S n +i is known 
exactly xa = 1/2. For chains up to length L = 204 and block Hilbert space size M — 256 we find results 

consistent with the theoretical prediction. Namely, if we calculate xa by fitting the log-log plot in Fig. || to a linear 
function, for points L > 10, we find xa ~ 0.6. The discrepancy with the exact result we attribute to the presence 
of a marginally irrelevant operator p6[ |, and other less relevant scaling operators. These corrections to scaling, which 
lead to a slightly larger value of the effective exponent, are evidenced by a slight curvature of the plot in Fig. ||. 
More precisely, the marginally irrelevant operator makes multiplicative logarithmic corrections to scaling p5| . We 
take these corrections into account by fitting the values to the functional form A(L/2) = L~ x& (log(L)) 7 /(L _1 ). 
The function fix) is assumed analytic around the origin and we expand it via Taylor series, keeping only the first 
two terms, f(x) = C± + C^x. With this procedure we obtain xa ~ 0.5 and 7 « 0.5, in agreement with theoretical 
expectations. 
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FIG. 2. Expectation value of the dimerization operator on the central bonds, A(L/2), for the S = 1/2 Heisenberg spin 
chain; L is the chain length. The infinite-size DMRG method with a block Hilbert space size of M — 256 is used. Slight 
curvature of the plot signals the presence of other less relevant scaling operators and logarithmic corrections to scaling. 

Returning to the supersymmetric chain, scaling is evident in Fig. [| where xa ~ 1-8 is calculated by fitting the 
M = 160 points with L > 10 to a straight line; curvature in the plot signals once again corrections to scaling. We 
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account for these less relevant operators by fitting A(L/2) to the functional form A(L/2) = L~ XA f(L~ 1 ). As in the 
case of the SU (2) Heisenberg chain above, we Taylor-expand f(x) = C\ + C2X + ■ ■ •, keeping only the first two terms. 
Following this procedure we obtain x& = 1-58 ± 0.07. The error is estimated by choosing different sets of points 
to include in the fit; statistical errors are smaller, at the few percent level. This value of the scaling dimension is 
consistent with the more precise value xa = 2 — v^ 1 = 1.57 ± 0.01 (using v — 2.33 ± 0.05), reported by other groups 
1- 

The DMRG works best for phases with a gap to low-lying excitations; we do not expect, nor do we find, very good 
quantitative agreement at the critical point. To obtain more accurate scaling it would be necessary to employ the 
"finite-size" DMRG algorithm and carry out careful finite-size analysis by varying both the chain length L and the 
dimerization parameter R to search for a universal scaling function. 
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FIG. 3. Expectation value of the dimerization operator on the central bonds, A(L/2), versus system length L for U — and 
at the critical point 7? = of the SUSY chain. Results for two different runs are shown: filled diamonds are for a block Hilbert 
space size of M = 160 [corresponding to a total Hilbert space size of (4M) 2 = 409, 600] and fi c = — 10 -6 in units where J = 1. 
Open circles are for a more accurate calculation at M = 256 [corresponding to a total Hilbert space size of 1,048,576] and 
He — — 10~ 4 . The two calculations agree precisely up through chain length L — 8. For L > 8 the entire Hilbert space cannot be 
accommodated; instead it is truncated via the DMRG algorithm. Nevertheless, the two calculations continue to agree closely 
out through chain length L — 22, the largest chain length reached in the M = 256 calculation. At L = 26 the M = 160 point 
deviates from the expected scaling form. We attribute the breakdown of scaling at large chain length to the truncation of the 
Hilbert space. 



VIII. DISCUSSION 



What have we learned from the above analysis? First, the scaling in Fig. || demonstrates the utility of the su- 
persymmetric approach to understanding transport at the IQHE plateau transition. The series of transformations 
outlined above which map the network model described by Eq. ^ to the pure, but rather complicated, supersymmetric 
spin-chain Hamiltonian Eq. ^5] appears to work. The replica trick is avoided and we obtain instead a concrete model 
amenable to analytical and numerical study. Deep inside a plateau (R — > ±1) the chain breaks up into disconnected 
two-site systems and Anderson localization in the interior is complete. At R = — 1 the two spins at the chain ends 
are fully polarized in the "[-direction and no extended edg e currents flow; in the opposite R = +1 limit the two spins 
are free and the conductivity increases by a quantum |27| ] . 

Second, and more importantly, the Hamiltonian is analytically tractable when projected onto the Hilbert space 
with only one type (f) of boson and fermion. This is the U 3> J limit, Eq. ^5], and there is the intriguing possibility 
that the Hamiltonian can be studied systematically down to the physical regime of finite U. This approach differs 
fundamentally from previous weak-coupling studies |^] where runaway RG flows toward strong coupling in the generic 
IQHE case suggested the existence of a strong-coupling fixed point. The starting point here, by contrast, would be at 
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very strong coupling. It may then be possible to incorporate spin excitations perturbatively about the strong coupling 
limit. 

A number of unanswered questions remain. The significance of the defectiveness of the Hamiltonian is unclear. 
It is a robust feature, as it occurs at both infinite and finite U. Also unanswered is the question of whether or not 
it will be possible to construct a two-dimensional conformal field theory (CFT) of the critical point. As pointed 
out by Mudry et al. p8f a candidate CFT would have to be non-unitary, thus allowing for scaling dimensions not 
predicted by the Kac determinant. This opens up the possibility of identifying a CFT that contains an operator with 
a scaling dimension xa ~ 11/7 in its spectrum; such an operator, if associated with a perturbation that induces a 
flow from the critical point to the non-critical phases of two adjacent Hall plateaus, would give rise to the expected 
correlation length exponent v = (2 — za)^ 1 ~ 7/3. Searching for such a field theory among CFT's which are also 
supersymmetric has not been successful p9[ . We believe that further studies of the supersymmetry approach to the 
Chalker-Coddington network model will shed more light on this matter. 
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